Thermoelectrical manipulation of nano-magnets 
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We investigate the interplay between the thermodynamic properties and spin-dependent transport 
in a mesoscopic device based on a magnetic multilayer (F/f/F), in which two strongly ferromag- 
netic layers (F) are exchange-coupled through a weakly ferromagnetic spacer (f) with the Curie 
temperature in the vicinity of room temperature. We show theoretically that the Joule heating pro- 
duced by the spin-dependent current allows a spin-thermo-electronic control of the ferromagnetic- 
to-paramagnetic (f/N) transition in the spacer and, thereby, of the relative orientation of the outer 
F-layers in the device (spin-thermo-electric manipulation of nanomagnets) . Supporting experimental 
evidence of such thermally controlled switching from parallel to antiparallel magnetization orien- 
tations in F/f(N)/F sandwiches is presented. Furthermore, we show theoretically that local Joule 
heating due to a high concentration of current in a magnetic point contact or a nanopillar can be used 
to reversibly drive the weakly ferromagnetic spacer through its Curie point and thereby exchange 
couple and decouple the two strongly ferromagnetic F-layers. For the devices designed to have an 
antiparallel ground state above the Curie point of the spacer, the associated spin-thermionic parallel- 
to-antiparallel switching causes magneto-resistance oscillations whose frequency can be controlled 
by proper biasing from essentially DC to GHz. We discuss in detail an experimental realization of 
a device that can operate as a thermo-magneto-resistive switch or oscillator. 



I. INTRODUCTION 

The problem of how to manipulate magnetic states 
on the nanometer scale is central to applied magneto- 
electronics. The torque effect 1,2 , which is based on the 
exchange interaction between spin-polarized electrons in- 
jected into a ferromagnet and its magnetization, is one of 
the key phenomena leading to current-induced magnetic 
switching. Current-induced precession and switching of 
the orientation of magnetic moments due to this effect 
have been observed in many experiments 3-12 . 

Current-induced switching is, however, limited by the 
necessity to work with high current densities. A natural 
solution to this problem is to use electrical point contacts 
(PCs). Here the current density is high only near the PC, 
where it can reach 13,14 values ~ 10 9 A/cm 2 . Since almost 
all the voltage drop occurs over the PC the characteristic 
energy transferred to the electronic system is comparable 
to the exchange energy in magnetic materials if the bias 
voltage Vbias ~ 0.1 V, which is easily reached in exper- 
iments. At the same time the energy transfer leads to 
local heating of the PC region, where the local tempera- 
ture can be accurately controlled by the bias voltage. 

Electrical manipulation of nanomagnetic conductors 
by such controlled Joule heating of a PC is a new prin- 
ciple for current-induced magnetic switching. In this pa- 
per we discuss one possible implementation of this prin- 
ciple by considering a thermoelectrical magnetic switch- 
ing effect. The effect is caused by a non- linear interac- 
tion between spin-dependent electron transport and the 
magnetic sub-system of the conductor due to the Joule 



heating effect. We predict that a magnetic PC with a 
particular design can provide both voltage-controlled fast 
switching and smooth changes of the magnetization di- 
rection in nanometer-size regions of the magnetic mate- 
rial. We also predict temporal oscillations of the magne- 
tization direction (accompanied by electrical oscillations) 
under an applied DC voltage. These phenomena are po- 
tentially useful for microelectronic applications such as 
memory devices and voltage controlled oscillators. 



II. EQUILIBRIUM MAGNETIZATION 
DISTRIBUTION 



The system under consideration has three ferromag- 
netic layers coupled to a non-magnetic conductor as 
sketched in Fig. 1. We assume that the Curie temper- 
ature T c (1) of region 1 is lower than the Curie temper- 
atures rj ' 2 ^ of regions and 2; in region 2 there is a 
magnetic field directed opposite to the magnetization of 
the region, which can be an external field, the fringing 
field from layer 0, or a combination of the two. We re- 
quire this magnetostatic field to be weak enough so that 
at low temperatures T the magnetization of layer 2 is 
kept parallel to the magnetization of layer due to the 
exchange interaction between them via region 1 (we as- 
sume the magnetization direction of layer to be fixed) . 
In the absence of an external field and if the tempera- 
ture is above the Curie point, T > Tc , the spacer of 
the proposed F/f(N)/F tri- layer is similar to the antipar- 
allel spin-flop 'free layers' widely used in memory device 
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FIG. 1: Orientation of the magnetic moments in a stack of 
three ferromagnetic layers (0, 1, 2) coupled to a non-magnetic 
conductor (3) ; the right arrow indicates the presence of a mag- 
netic field H , which is antiparallel to the stack magnetization. 




FIG. 2: Sketch of the spatial dependence of the orientation of 
the magnetic moments in the stack of Fig. 1 at a temperature 
T higher than the temperature Tj° r ', at which the parallel 
orientation becomes unstable, but lower than the Curie tem- 
perature Tj 1 ^ of layer (1). 



applications 1 '. 

As T approaches Tc from below the magnetic mo- 
ment of layer 1 decreases and the exchange coupling be- 
tween layers and 2 weakens. This results in an inhomo- 
geneous distribution of the stack magnetization, where 
the distribution that minimizes the free energy of the 
system is given by Euler's equation (see, e.g., Ref. 16): 



4- (a{ X )M\ X )f 
ax V ax 



■^-M 2 sin 28+— S m9 



Here the x-axis is perpendicular to the layer planes of the 
stack, the z-axis is directed along the magnetization di- 
rection in region 0; the magnetization direction depends 
only on x and the vector M rotates in-plane (that is in 
the yz-plane) 16 ; 9{x) is the angle between the magnetic 
moment M(x) at point x and the z-axis (in the yz-planc) 
and M{x) = |M(x)|. In the case under consideration 
<x{x) = ax, (3{x) = f3\ for < x < L\ and a(x) = a 2 , 
P{x) = /3 2 for L\ < x < L 2 ; here ai j2 ~ /i^/aM^, 

(1 2) 

where a is the lattice spacing, 

h 2 ~ ksTc and Mi 2 
are the exchange energies and magnetic moments of re- 
gions 1 and 2; /3i (/3 2 ) is a dimensionless measure of the 
anisotropy energy of region 1 (region 2); ks is the Boltz- 
mann constant. Below we assume the lengths Li j2 of re- 
gions 1 and 2 to be short er than th e domain wall lengths 
in these regions Z1.2 = -\/c*i,2//3i, 2 . 

In order to find the magnetization distribution inside 
the stack one may solve Eq. (1) in regions 1 and 2 to 



get 8i(x) and 8 2 (x), respectively, and then match these 
solutions at the magnetization interface x = L\. Inte- 
grating Eq. (1) with respect to x in the limits L\ — 6 < 
x < L\ + S, 5—^0 one gets the matching condition as 
follows: 



a 2 M, 



,d9 2 (x) 



dx 



x—Li 



aiM( 



2 d6i(x) 



d/x 



x—Li 



9 2 (L 1 ) = 9 1 (Li). (2) 

The boundary condition at the ferromagnetic interface 
x = between layers and 1 follows from the require- 
ment that the direction of the magnetization in layer is 
fixed along the z-axis (i.e., 0(x) = in this layer): 



0i(O) =0. 



(3) 



At the "free" end of the ferromagnetic sample the 
boundary condition for the magnetization M(x) is 
dM{x)/dx = (see, e.g. Ref. 17), so that 



d6 2 {x) 
dx 



0. 



(4) 



Solving Eq. (1) in regions 1 and 2 under the assump- 
tion L lj2 <g; Z lj2 and with the boundary conditions (2) - 
(4) one finds the magnetization in region 1 to be inho- 
mogencous, 



Li «i 



< x < Li, 



(5) 



while due to the boundary condition (4) the magnetic 
moments in region 2 are approximately parallel, to within 
corrections of order aiM 1 2 (T)L 2 /a 2 A/|(T)Li <C 1, i.e. 



0- (1) 2 (x) = e 2 (L 1 +L 2 ) 



sin 6», (Li 



8a 2 M 2 



■(L. + L2-X) 2 



(6) 

where L± < x < L\ + L 2 . Using the above boundary 



conditions one finds that Q 2 (L\) 
determined by the equation 

6 = D(H, T) sin6, 

e = ±tt, 



l {L l + L 2 ) = 9 is 



T < T c (1) 

t > ri 1 ' 



where 



D(H,T) 



LiL 2 HM 2 (T) 
4a 1 M 1 2 (T) 



(7) 



(8) 



In Eq. (8) Mi(T) = M?>^J{TP -T)/T^> and M 2 (T) 
are the magnetic moments of region 1 and 2, respectively; 
the parameter D(H, T) is the ratio between the magnetic 
energy and the energy of the stack volume for the in- 
homogeneous distribution of the magnetization. As the 
second term inside the brackets in Eq. (6) is negligibly 
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small, the magnetization tilt angle 9 in region 2 becomes 
independent of position and is simply given as a function 
of H and T by Eq. (7). 

By inspection of Eq. (7) one finds that it has either one 
or several roots in the interval — n < O < 7r depending 
on the value of the parameter D(H, T). 

At low temperatures the exchange energy prevails, 
the parameter D(H,T) < 1 and Eq. (7) has only one 
root, = 0. Hence a parallel orientation of all mag- 
netic moments in the stack is thermodynamically sta- 
ble. However, at temperature Tj 01 - 1 < Tc for which 
D(Tc , H) = 1, two new roots = ±|# m i n | ^= appear. 
The parallel magnetization corresponding to = is 
now unstable 18 and the direction of the magnetization in 
region 2 tilts as indicated in Fig. 2. Using Eq.(8) one finds 
the critical temperature of this orientation transition to 
be equal to 



T (or) = T (i) 1 



5T L X L 2 HM 2 



ST 



T^J' T c (1) 4a 1 M 1 2 (0) 



(9) 

The tilt increases further with T until at T = T c (1) the 
exchange coupling between layers and 2 vanishes and 
their magnetic moments become antiparallel. 



Thermally assisted exchange decoupling in F/f/F 
multilayers 



ature leads at the same time to a lower switching field of 
the magnetically hard layer, which is due to the stronger 
effective magnetic torque on the top layer in the coupled 
F0/f/F2 state. This thermally-controlled interlayer ex- 
change coupling is perfectly reversible on thermal cycling 
within the given temperature range. 
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FIG. 3: Magnetization loop for a sample of SiC>2/Cu 
90/Ni 80 Fe 2 o 8/Co go Fe 10 2/Ni 7 oCu 3 o 30/Co 90 Feio 5/Ta 10 
[nm] as the temperature is varied from 25° C to 130° C. Hco 
and Hc2 are the coercive fields of the bottom and top mag- 
netic layers, respectively. 



To demonstrate the properties of the tri-layer mate- 
rial system proposed above we have suitably alloyed Ni 
and Cu to obtain a spacer with a T c just higher than 
room temperature (RT). The alloying was done by co- 
sputtering Ni and Cu at room temperature and base pres- 
sure 10 -8 torr on to a 90x10 mm long Si substrate in such 
a way as to obtain a variation in the concentration of Ni 
and Cu along the substrate. By cutting the substrate 
into smaller samples along the compositional gradient, a 
series of samples were obtained, each having a different 
Curie temperature. 

One of the multi-layer compositions chosen was NiFe 
8/CoFe 2/NiCu 30/ CoFe 5 [nm] where the NiFe layer is 
used to lower the coercive field of the bottom layer (Hco) 
in order to separate it from the switching of the top layer 
(Hc2)- A magnetometer equipped with a sample heater 
was used to measure the magnetization loop as the tem- 
perature was varied between 25°C and 130°C. The re- 
sults for a Ni concentration of ^70% are shown in Fig. 3. 
The strongly ferromagnetic outer layers are essentially 
exchange- decoupled at T > 100°C (F/Paramagnetic/F 
state), as evidenced by the two distinct magnetization 
transitions at approximately 15 and 45 Oe in Fig. 3. As 
the temperature is reduced to RT, the switching field of 
the soft layer increases and the originally sharp M-H 
transition becomes significantly skewed. This confirms 
the theoretical result, expressed by Eqs. (5) and (6) for 
Q(H,T), that the magnetic state of the sandwich is of 
the spring-fcrromagnct type 1 ' 1 . The lowering of tcmpcr- 



We further demonstrate an exchange-biased magnetic 
tri-layer of the generic composition AF/F0/f/F2, where 
the spacer separating the outer ferromagnetic layers (F) 
is a low-Curie temperature diluted ferromagnetic alloy (f ) 
and one of the F0 layers is exchange-pinned by an antifcr- 
romagnet (AF). In addition to the tri-layer a Cu spacer 
and a reference layer, pinned by an AF, have been added 
on top of the stack in order to measure the current-in- 
plane giant magnetoresistance (GMR). The specific stack 
composition chosen was Si/Si02/NiFe 3/MnIr 15/CoFe 
2/ Ni 70 Cu 30 30/CoFe 2/NiFe 10/CoFe 2 /Cu 7/CoFe 
4/NiFc 3/MnIr 15/Ta 5 [nm]. The sample was deposited 
at room temperature in a magnetic field of 350 Oe, then 
annealed at 300°C for 20 minutes, and field cooled to 
RT in ~ 800 Oe. The NiCu spacer was co-sputtered 
while rotating the substrate holder, such that the final 
concentration was 70% Ni and 30% Cu having the Tq 
suitably above RT. Fig. 4 shows how the interlayer 
exchange field H ox of this sample varies with tempera- 
ture. H ea; shown in the main panel of Fig. 4 is de- 
fined as the mid point switching field of the soft F2- 
laycr (~ 18, - 32, and ~ 47 Oe for 100°C, 60°C, and 
25°C, respectively; see inset), which reflects the strength 
of the interlayer exchange coupling through the spacer 
undergoing a ferromagnetic-paramagnetic transition in 
this temperature range. To explain why this is so, we 
need to consider the difference in effective magnetic thick- 
ness between the top and bottom pinned ferromagnetic 
layers. The effective magnetic thickness for the bot- 



4 



torn pinned CoFe/NiCu/CoFe/NiFe/CoFe layers is ap- 
proximately three times larger than for the top pinned 
CoFe/NiFe. From the inset to Fig 4, the temperature 
variation of the exchange pinning for the top pinned 
CoFe/NiFe is 20 Oe or 0.3 Oe/K. If we were to assume 
that the bottom pinned CoFe/NiCu/CoFe/NiFe/CoFe 
layers are coupled and reverse as one layer, and that the 
variation in exchange field is caused solely by the weaken- 
ing pinning at the bottom Mnlr interface, then we would 
expect an exchange field three times smaller than for the 
top pinned CoFe/NiFe. With a three times smaller ex- 
change field at RT the expected temperature variation 
would be 7 Oe or 0.1 Oe/K, which clearly is much lower 
than the observed change of 25 Oe (from 45 Oe to 20 Oe) 
and therefore the measured de-pinning of the switching 
layer is predominantly due to a softening of the exchange 
spring. 

We have separately measured the strength of the ex- 
change pinning at the bottom Mnlr surface. For CoFe 
ferromagnetic layers 2-4 nm thick, the pinning strength 
at RT is 500 Oe or more. At 130 C, at which the spacer 
is paramagnetic and fully decoupled from the underlying 
Mnlr/CoFe bilayer, the pinning strength is still above 100 
Oe. We therefore conclude that the dominating effect in 
question is the weakening exchange spring in the spacer. 
This demonstrates the principle of the thermionic spin- 
valve proposed, where the P to AP switching is controlled 
by temperature. The AF-pinncd implementation of the 
spin-thermionic valve presented should be highly relevant 
for application. 
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FIG. 4: Interlayer exchange field H ex versus temper- 
ature T. The composition of the complete stack is 
Si/SiOa/NigoFeao 3/ Mn 80 Ir 20 15/Co 9 oFe 10 2/ Ni 70 Cu 3 o 
30/Co 90 Fe 10 2/Ni 80 Fe 2 olO/Co 9 oFe 10 2 /Cu 7/Co 90 Fe 10 
4/NigoFe2o 3/Mn8oIr20 15/Ta 5 [nm]. Inset: Current-in-plane 
GMR at T=25, 60 and 100° C. 

As is obvious from the above analysis the dependence 
of the magnetization direction on temperature allows 
electrical manipulations of it by Joule heating with an 



applied current flowing through the stack. In the next 
section we find connection between the magnetization di- 
rection and the current-voltage characteristics (IVC) of 
such a spin-thermionic valve. 



III. THERMOELECTRIC MANIPULATION OF 
THE MAGNETIZATION DIRECTION. 

A. Current-voltage characteristics of the stack 
under Joule heating. 

If the stack is Joule heated by a current J its temper- 
ature T(V) is determined by the heat-balance condition 

JV = Q{T) 1 J = V/R{&), (10) 

and Eq. (7), which determines the temperature depen- 
dence of Q(T(V)). Here Q(T) is the heat flux from the 
stack and R(Q) is the stack resistance. In the vicinity of 

the Curie temperature T c Eq. (7) can be re-written as 

f ±tt, T > T c (1) 

9= L rw . T ~ (i) (11) 

[ A) T (i)_ T sin0, T < Tc , 

(here D is defined in Eq.(9)). 

Equations (10) and (11) define the current-voltage 
characteristics (IVC) of the stack, J = G(Q(V))V, G = 
R , in a parametric form which can be re- written as 

V = vVc (1) )^)(l-^^) . (12) 

The parameter 9 is defined in the interval —tt<9<tt, 
T dQ, 

D = DO Qdlf\T=TP^ D 0, 

and in order to derive Eq. (12) we used the expansion 
Q(T) = Q(T C (1) ) + Q' T {T^){T - T c (1) ) [Q' T = dQ/dT}. 

It follows from Section II that the stack resistance is 
R(0) in the entire temperature range T(V) < Tc and 

R(n) in the range T(V) > T c (1) . This implies that the 
IVC branches J — G(Q)V and J = G(ir)V are linear for, 
respectively, 

V<Vi = \j i?(0)(5(T c (or) ) (13) 
(0 — a in Fig. 5) and 

V>V C = \/ R{tt)Q(T^ 1} ) (14) 

(b — b' in Fig. 5). If V% < V < V c the stack temperature 
is T c (or) < T(V) < T c (1) , and the direction of the magne- 
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tization in region 2 changes with a change of V; hence 
the IVC is non-linear there. Below we find the condi- 
tions under which this branch of the IVC has a negative 
differential conductance. 

Differentiating Eq. (12) with respect to V one finds 



dJ ffffl) [G(fl)(l 



D sin 6/6)]' 



D sin 0/0)]' 



e=e(v) 



(15) 



where [...]' means the derivative of the bracketed quan- 
tity with respect to the angle 0, and Q(V) is found from 
the second equation in Eq. (12). From this result it fol- 
lows that the differential conductance Gd{V) = dJ/dV is 
negative if 



d (1 -DsinO/e) 
d& i?(9) 

For a stack resistance of the form 



< 0. 



where 



R(Q) = R+(l-rcose) 

J?_ RM ± R(0) 
r= iV R± = 2 



(16) 



(17) 



one finds that the differential conductance dJ/dV < if 

3r 



D < 



1 + 2r 



(18) 



Hence the IVC of the stack is N-shaped as shown in 
Fig. 5. 

We note here that the modulus of the negative differ- 
ential conductance may be large even in the case that the 
magnetoresistance is small. Using Eq.(15) at r < 1 one 
finds the differential conductance Gdiff as 



G 



diff 



dJ 

~dV 



-R-\0) 



1 - Do/'Sr 
1 + D Q /3r 



(19) 



which is negative provided Dq < 3r, the modulus of 
Gdiff being of the order of i? _1 (0). 

Here and below we consider the case that the elec- 
tric current flowing through the sample is lower than 
the torque critical current and hence the torque effect 
is absent 20 

As the IVC curve J(V) is N-shaped the thermoelec- 
trical manipulation of the relative orientation of layers 
and 2 may be of two different types depending on the 
ratio between the resistance of the stack and resistance 
of the circuit in which it is incorporated 

In the voltage-bias regime which corresponds to the 
case that the resistance of the stack is much larger than 
resistance of the rest of the circuit, the voltage drop 
across the stack preserves the given value which is ap- 
proximately equal to the bias voltage and hence there 
is only one value of the current (one point on the IVC) 
J = Juas corresponding to the bias voltage Vuas (see 
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FIG. 5: Current-voltage characteristics (IVC) of the mag- 
netic stack of Fig. 1 calculated for R(0) = R+ — R- cos O, 
R-/R+ = 0.2, D = 0.2; J c = V c /R{tv). The branches - a 
and b — b' of the IVC correspond to parallel and antiparal- 
lel orientations of the stack magnetization, respectively (the 
parts a — a' and — b are unstable); the branch a — b cor- 
responds to the inhomogeneous magnetization distribution 
shown in Fig. 2. 
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FIG. 6: The angle 0, which describes the tilt of the direction 
of the magnetization in layer 2 with respect to that in layer 
1 (see Fig. 2), as a function of voltage in the voltage- biased 
regime (top) and current in the current-biased regime (bot- 
tom). Both curves were calculated for R(Q) = R+—R- cos Q, 
R-/R+ = 0.2, D = 0.2; J c = V c /R{tt). 
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Fig. 5) In this case the relative orientation of the magne- 
tization of layers and 2 can be changed smoothly from 
being parallel to anti-parallel by varying the bias voltage 
through the interval V\ < Vbias _ V c . This corresponds 
to moving along the a — b branch of the IVC. The depen- 
dence of the magnetization direction on the voltage 
drop across the stack is shown in Fig. 6. 

In the current-bias regime, on the other hand, which 
corresponds to the case that the resistance of the stack 
is much smaller than the resistance of the circuit, the 
current in the circuit J is kept at a given value which 
is mainly determined by the bias voltage and the circuit 
resistance (being nearly independent of the stack resis- 
tance). As this takes place, the voltage drop across the 
stack V differs from the bias voltage Vbias, being deter- 
mined by the equation J(V) = J. As the IVC is re- 
shaped, the stack may now be in a bistable state: if the 
current is between points a and b there are three possible 
values of the voltage drop across the stack at one fixed 
value of the current (see Fig. 5). The states of the stack 
with the lowest and the highest voltages across it are sta- 
ble while the state of the stack with the middle value of 
the voltage drop is unstable. Therefore, a change of the 
current results in a hysteresis loop as shown in Fig. 6: 
an increase of the current along the — a' branch of the 
IVC leaves the magnetization directions in the stack par- 
allel (0 = 0) up to point a, where the voltage drop V 
across the stack jumps to the right branch b — b', the 
jump being accompanied by a fast switching of the stack 
magnetization from the parallel to the antiparallel orien- 
tation (0 = ±7r). A decrease of the current along the 
6' — IVC branch keeps the stack magnetization antipar- 
allel up to point 6, where the voltage jumps to the left 
— a' branch of the IVC and the magnetization of the 
stack comes back to the parallel orientation (0 = 0). 

In the next Section we will show that this scenario for 
a thermal-electrical manipulation of the magnetization 
direction is valid for small values of the inductance in the 
electrical circuit. If the inductance exceeds some critical 
value the above steady state solution becomes unstable 
and spontaneous oscillations appear in the values of the 
current, voltage drop across the stack, temperature, and 
direction of the magnetization. 



B. Self-excited electrical, thermal and directional 
magnetic oscillations. 

1. Current perpendicular to layer planes (CPP) 

Consider now a situation in the bias voltage regime 
where the magnetic stack under investigation is con- 
nected in series with an inductance C and biased by a 
DC voltage Vbi as , as described by the equivalent circuit 
in Fig. 7. The thermal and electrical processes in this 




FIG. 7: Equivalent circuit for a Joule-heated magnetic stack 
of the type shown in Fig. 1. A resistance R(V) = J(t)/V(t), 
biased by a fixed DC voltage Vbias, is connected in series with 
an inductance £; V(t) is the voltage drop over the stack and 
J(t) is the total current. 



system are governed by the set of equations 

C v %- = J 2 R(e)-Q(T); C^- + JR(@) = V hias , (20) 
at at 

where Cy is the heat capacity. The relaxation of the 
magnetic moment to its thermodynamically equilibrium 
direction is assumed to be the fastest process in the prob- 
lem, which implies that the magnetization direction cor- 
responds to the equilibrium state of the stack at the 
given temperature T(t). In other words, the tilt angle, 
= Q(T(t)), adiabatically follows the time-evolution of 
the temperature and hence its temperature dependence 
is given by Eq. (7). 

A time dependent variation of the temperature is ac- 
companied by a variation of the magnetization angle 
Q(T(t)) and hence by a change in the voltage drop across 
the stack via the dependence of the magneto-resistance 
on this angle, R = R(Q). 

The system of equations Eq.(20) has one time- 
independent solution (T '(Vbias), J (Vbias)) which is deter- 
mined by the equations 

J 2 R (0(T)) = Q(T), JR(&(T)) = Vbias (21) 

This solution is identical to the solution of Eqs.(7,10) 
that determines the N-shaped IVC shown in Fig. 5 with 
a change J — > J and V — > Vbi as ■ 

In order to investigate the stability of this time- 
independent solution we write the temperature, current 
and the angle as a sum of two terms, 

T = f (Vuas) + Ti(t); 
J = J (Vbias) + Ji(t); 

Q = ®(Vbias) + 0l(t), (22) 

where T\ , Ji and 9\ each is a small correction. Insert- 
ing Eq.(22) into Eq.(20) and Eq. (7) one easily finds that 
the time-independent solution Eq.(21) is always stable 
at any value of the inductance C if the bias voltage Vbias 
corresponds to a branch of the IVC with a positive differ- 
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inductance exceeds the critical value Eq.(24) the system 
looses its stability and a limit cycle appears in the plane 
(J, T) (see, e.g., 22 ). This corresponds to the appearance 
of self-excited, non-linear and periodic temporal oscil- 
lations of the temperature T = T(t) and the current 
J = J(t), which are accompanied by oscillations of the 
voltage drop across the stack V(t) = J(t) and the the 
magnetization direction 0(f) = 0(T(t)). For the case 
that (£ — C c )/C c -C 1 the system executes nearly har- 
monic oscillations around the steady state (see Eq.(22)) 
with the frequency u> = Im7(£ = C c ), that is the tem- 
perature T, the current J, the magnetization direction 
and the voltage drop across the stack V(t) = R(Q(t))J(t) 
execute a periodic motion with the frequency 



(25) 



With a further increase of the inductance the size of 
the limit cycle grows, the amplitude of the oscillations 
increases and the oscillations become anharmonic, the 
period of the oscillations therewith decreases with an in- 
crease of the inductance C. 

In order to investigate the time evolution of the voltage 
drop across the stack and the current in more details 
it is convenient to introduce an auxiliary voltage drop 
V(t) and a current Jo(t) related to each other through 
Eqs. (10) and (11). Hence we define 



FIG. 8: Spontaneous oscillations of the current J(t) and the 
voltage drop V(t) over the stack calculated for R-/R+ — 0.2, 
Do = 0.2 and (£ - £ C r)/£cr = 0.013; J c = V c /R(tv). J(t) and 
V(t) develop from the initial state towards the limit cycle 
(thick solid line) along which they execute a periodic motion. 
The thin line is the stationary IVC of the stack. The bottom 
figure shows the limit cycle along which Q(t) and V(t) execute 
a periodic motion. 



cntial resistance (branches 0-a and b-b' in Fig. 5). If the 
bias voltage Vu as corresponds to the branch with a neg- 
ative differential resistance (Vi < Vuas < V c , see Fig. 5) 
the solution of the set of linearized equations is T\ = 



t[ 0) cxp{ 7 t}, Ji = J x w cxp{ 7 t} and 6 



(0), 



where T±°\ J^' and 9 { "' are any initial values close to 
the steady-state of the system, and 



7(0) 



?i 0) exp{ 7 i} 



R /C-C c 
2£V C r 



± 



where 



Cc = 



'C-C c 



CV 



\d(GQ)/dT\ 



T=T(V) 



(24) 



and Rd = dV/dJ, R = i?(0) is the differential resistance. 

As is seen from Eq.(23) the steady-state solution 
Eq.(21) is stable only if the inductance C < C c ; if the 



V(t) = y/R(T(t))Q(T(t))); J = V(t)/R(T(t)) , (26) 

where R(T) = R(Q(T)). Comparing these expressions 
with Eq. (10) one sees that at any moment t Eq. (26) 
gives the stationary IVC of the stack, J = Jq{V), defined 
by Eq. (12) (changing J — >• J and V -> V), see Fig. (5). 

Differentiating V(t) with respect to t and using 
Eqs. (20) and (26) one finds that the dynamical evolution 
of the system is governed by the equations 

dV J 2 - 4{V) 



T0 H = 



2J (V) 



„d.J JV 
C— H =r- 

dt J (V) 



(27) 



where 



To 



C v 



(QR)' T 



T=T(V) 



As follows from the second equation in Eq. (26), at 
any moment t the voltage drop over the stack V(t) = 
R(T(t))J(t) is coupled with V(t) by the following rela- 
tion: 



V 



J 



The coupled equations (20) have only one steady-state 
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FIG. 9: Spontaneous oscillations of the current J(t), the volt- 
age drop V(t), the magnetization direction angle O(t), and the 
temperature T(t) corresponding to motion along the limit cy- 
cle shown in Fig. 8. Calculation parameters are R-/R+ — 0.2, 
D = 0.2 and (£ - £ cr )/£cr = 0.3 x 10~ 4 ; J c = V c /R{tt). 



solution J = Jo(Vbias) where Jq(V) is the IVC shown 
in Fig. 5 (see Eqs. (26). However, in the interval V\ < 
Vbms < V c this solution is unstable with respect to small 
perturbations if £ > £ cr . As a result periodic oscilla- 
tions of the current J(t) and V(t) appear spontaneously, 
with J(t), V(t) eventually reaching a limit cycle. The 
limiting cycle in the J-V plane is shown in Fig. 8. The 
stack temperature T = T(t), the magnetization direction 
8(t) = 0(T(i)), follow these electrical oscillations adia- 
batically according to the relations Q(T(t)) = V(t)J (t) 
(here J (t) = J (V(t))) and B(t) = 8(T(i)) (see Eq. (7)) 
as shown in Fig. 9. 

The character of the oscillations changes drastically in 
the limit £ 3> £ C r- In this case the current and the volt- 
age slowly move along the branches — a and b — b of 
the IVC at the rate J/ J w R+/C, quickly switching be- 
tween these branches at the points a and b with the rate 
~ 1/to (see Fig. 10). Therefore, in this case the stack pe- 
riodically switches between the parallel and antiparallcl 
magnetic states (see Fig. 11). 



2. Current in the layer planes (CIP). 

If the electric current flows in the plane of the lay- 
ers (CIP) of the stack the torque effect is insufficient or 
absent 1 '- 1 while the magneto-thermal-electric oscillations 
under consideration may take place. In this case the to- 
tal current flowing through the cross-section of the layers 
may be presented as 

Jcip = {R-\e) + Ro 1 )V (28) 

where i?(0) and Rq are the magneto-resistance and the 
angle-independent resistance of the stack in the CIP set 
of the experiment. 

In a CIP configuration the stack is Joule heated by 
both the angle-dependent and the angle-independent cur- 




FIG. 10: Spontaneous oscillations of the current J(t) and the 
voltage drop V(t) calculated for R-/ R+ = 0.2, Do = 0.2 and 
(C — £ C r)/£cr = 535; J c = V c /R(n). The time development of 
J(t) and V(t) follows one or the other of the dashed lines to- 
wards the limit cycle (thick solid line) depending on whether 
the initial state is inside or outside the limit cycle. The bot- 
tom figure shows how the current oscillations develop if the 
initial state is inside the limit cycle. The stationary IVC of 
the stack is shown as a thin solid line. 
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FIG. 11: Spontaneous oscillations of the magnetization direc- 
tion angle 0(t) calculated for R-/R+ — 0.2, Do = 0.2 and 
(£ - £ C r)/£or = 535; J c = V c /R(tt). 
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rents and hence Eq.(lO) should be re- written as follows: 
JcipV = Q(T), J = V/R eff (e), (29) 

where 

R(Q)R 



i?(6) + Rq 



(30) 



Using Eq.(15) and Eq.(28) one finds that the pres- 
ence of the angle- independent current in the stack modi- 
fies the condition of the negative differential conductance 
dJcip/dV: it is negative if 



D < 



3r 



R -{l-rfR A 



(31) 

As is seen from here, an IVC with a negative differential 
resistance is possible if R a > (1 — r) 2 R + (sec Eq.[17] for 
definitions of R± and r) . 

The time evolution of the system is described by the 
set of equations Eq.(20) in which one needs to change 
J — > Jcip and i?(6) — > R e ff(Q). Therefore, under 
this change, the temporal evolution of the system in 
a CIP configuration is the same as when the current 
flows perpendicular to the stack layers: if the bias volt- 
age corresponds to the negative differential conductance 
dJcip/dV < and the inductance exceeds the critical 
value 



\d(G eff Q)/dT\ 



T=T(V) 



(32) 



where G e f / = R^ff ! self-excited oscillations of the cur- 
rent J civ, voltage drop over the stack V, the tempera- 
ture T and the angle 8(T(U)) arise in the system, the 
maximal frequency of which being 



y/\dV/dJ C Ip\Reff(T(V)) 



v=v bia 



(33) 



if (C - C cr )/C CI « 1 

Below we present estimations of the critical inductance 
and the oscillation frequency which are valid for both 
the above mentioned CPP and CIP configurations of the 
experiment. 

Using equations Eq.(21) and Eq.(24) one may esti- 
mate the order of magnitude of the critical inductance 
and the oscillation frequency as C c « Tc v /j 2 d and 
uj ~ pj 2 /Tc v where c v is the heat capacity per unit vol- 
ume, p is the resistivity, and d is a characteristic size of 
the stack. For point contact devices with typical values 
of d ~ 1CT 6 lCT 5 cm, c v ~ 1 J/cm 3 K, p ~ 1CT 5 Qcm, 
j ~ 10 s A/cm 2 and assuming that cooling of the device 

can provide the sample temperature T w ~ 10 2 K 
one finds the characteristic values of the critical induc- 
tance and the oscillation frequency as C cr ~ 10 _8 -^10~ 7 H 
and to » 1GHz 



IV. CONCLUSIONS. 

The experimental implementation of the new princi- 
ple proposed in this paper for the electrical manipula- 
tion of nanomagnetic conductors by means of a controlled 
Joule heating of a point contact appears to be quite fea- 
sible. This conclusion is supported both by theoretical 
considerations and preliminary experimental results, as 
discussed in the main body of the paper. Hence we ex- 
pect the new spin-thermo-electronic oscillators that we 
propose to be realizable in the laboratory. We envision 
F0/f/F2 valves where two strongly ferromagnetic regions 
(T c ~ 1000 K) are connected through a weakly ferromag- 
netic spacer (T c <C 1000 K). The Curie temperature of 
the spacer would be variable on the scale of room temper- 
ature, chosen during fabrication to optimize the device 
performance. For example, doping Ni-Fe with ~ 10% of 
Mo brings the T c from - 1000 K to 300-400 K. Alter- 
natively, alloying Ni with Cu yields a spacer with a T c 
just above Room Temperature (at RT or below RT, if 
needed). If a sufficient current density is created in the 
nano-tri-layer to raise the temperature to just above the 
T c of the spacer, the magnetic subsystem undergoes a 
transition from the F0/f/F2 state to an F0/N/F2 state, 
the latter being similar to conventional spin-valves (N for 
nonmagnetic, paramagnetic in this case). Such a transi- 
tion should result in a large resistance change, of the same 
magnitude as the "giant magnetoresistance" (GMR) for 
the particular material composition of the valve. 

Local heating (up to 1000 K over 10-50 nm) can read- 
ily be produced using, e.g., point contacts in the thermal 
regime, with very modest global circuit currents and es- 
sentially no global heating 1 '. Heat is known to propagate 
through nm-sized objects on the ns time scale, which can 
be scaled with size to the sub-ns regime. When voltage- 
biased to generate a temperature near T c (f), such a 
F0/f/F2 device would oscillate between the two magnetic 
states, resulting in current oscillations of a frequency that 
can be tuned by means of connecting a variable induc- 
tance in series with the device. Spin rotation frequencies 
may be tuned from the GHz-range down to quasi-DC (or 
DC as soon as the inductance is smaller than the critical 
value). For F0/f/F2 structures geometrically designed 
in the style of the spin-flop free layer of today's magne- 
toresistive random access memory (MRAM) , the dipolar 
coupling between the two strongly ferromagnetic layers 
would make the anti-parallel state (F0 t /N/F2 \. ) the 
magnetic ground state above T c (f). The thermal transi- 
tion in the f- layer would then drive a full 180-degree spin- 
flop of the valve. The proposed spin-thermo-electronic 
valve can be implemented in CPP as well as CIP geome- 
try, which should make it possible to achieve MR signals 
of 10. 

In conclusion, we have shown that Joule heating of the 
magnetic stack sketched in Fig. 1 allows the relative ori- 
entation of the magnetization of the two ferromagnetic 
layers and 2 to be electrically manipulated. Based 
on this principle, we have proposed a novel spin-thermo- 
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electronic oscillator concept and discussed how it can be 
implemented experimentally. 
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